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Abstract. We summarize our model for long-duration gamma ray bursts (GRBs) that fits the redshift (z) distributions mea- 
sured with Swift and missions before Swift, and the pre-Swift GRB jet opening-angle distribution inferred from achromatic 
breaks in the optical light curves. We find that the comoving rate density of GRB sources exhibits positive evolution to z ^ 3 
- 5, whereas the star formation rate inferred from measurements of the blue and UV luminosity density peaks at z ~ 1 - 3. 
The mean intrinsic beaming factor of GRBs is found to be « 34 - 42, and we predict that the mean GRB optical jet opening 
half-angle measured with Swift is « 10°. We estimate the number of GRBs per year that GLAST is expected to observe based 
on ratios of BATSE and EGRET GRB fluences. 
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> . INTRODUCTION 

(N : 

With the launch of the Swift satellite 1 1 ], rapid follow-up studies of GRBs triggered by the Burst Alert Telescope (BAT) 
on Swift became possible. A fainter and more distant population of GRBs than found with the pre-Swift satellites 
CGRO-BATSE, BeppoSAX, INTEGRAL, and HETE-2 is detected [2]. The mean redshift of 41 pre-Swift GRBs that 
also have measured optical afterglow beaming breaks [J] is (z) = 1.5, while 16 GRBs discovered by Swift have 
(z) = 2.72 I4|]. The null hypothesis that the redshift distributions of Swift and pre-Swift GRBs are the same is rejected 
in several statistical tests at the ^97% confidence level [5D. 

In recent work [6], we considered whether the differences between the pre-Swift and Swift redshift distributions can 
be explained with a physical model for GRBs that takes into account the different flux thresholds of GRB detectors, 



o 
t> 
o 

Oh; 

where we assume the detector threshold for Swift and pre-Swift GRBs to be ~ 10 -8 and ~ 10 -7 ergs cm -2 s~ 
+i . respectively [6], corresponding to the minimum peak fluxes of GRBs with redshifts measured by Swift and missions 

before Swift. We summarize J60 and make an independent estimate of the number of GRBs that GLAST should detect. 
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: RESULTS 

We parameterize the distribution of GRB jet opening angles with a function of the form dN/djij oc (1 — fij)\ where 
arccos(/iy) is the jet half-angle, and find best fit values for the y-ray energy release §y for different functional forms 
of the comoving rate density of GRBs (see Fig. 1, left panel). We simplify the analysis with a flat GRB vF v spectrum, 
and assume that the properties of GRBs do not change with time. Adopting the uniform jet model where the energy 
per solid angle is roughly constant throughout the GRB jet, and taking a mean intrinsic duration of 10 seconds, we 
obtain best-fit values with a corrected emitted y-ray energy (f r = 4x 10 51 ergs, jet opening angles in the range 0.05 - 
0.7 radians, and with s « —1.25. 

Our analysis shows that good fits to the pre-Swift and Swift redshift and opening angle data require a GRB rate 
history that rises faster than the star formation rate (SFR2,3,4) at high redshifts (SFR5 and SFR6; see Fig. 1, left 
panel). The results of our fitting indicate that GRB activity was greater in the past and is not simply proportional 
to the bulk of the star formation as traced by the blue and UV luminosity density of the universe. Furthermore, our 
model predicts the mean intrinsic beaming factor of GRB jetted outflows contributing to the optical breaks to be in 
the range from « 34 — 42. The analysis also indicates that the average jet opening half-angle of GRBs detected with 
Swift is (0j) ~ 10°, compared to the pre-Swift average of 7° |6]. Thus we expect to detect more faint low-redshift, 
large opening angle GRBs that pre-Swift satellites could not detect (see also Ji). 

We also estimate the number of GRBs per year GLAST will observe using the fluence ratios of BATSE (~ 20 keV 
- 2 MeV) and EGRET (100 MeV - 5 GeV) for the 5 BATSE bursts also detected in the EGRET spark chamber U. 
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FIGURE 1. Left panel: star formation rate (SFR) functions used in our GRB study. The solid line (SFR1) is a constant comoving 
density; SFR2 and SFR4 are low and high ranges for the SFR; SFR3 is a recent fit 0] to the data; SFR5 and SFR6 are GRB rate 
histories that give a good fit to the Swift and pre-Swift redshift distribution and the pre-Swift opening angle distribution |6]. Right 
panel: Estimate of the number of GRBs per year GLAST will observe, based on measured EGRET/BATSE fluence ratios. 



We find that the average BATSE/EGRET fluence ratio is « 5%, though with large scatter and not taking into account 
Earth occultation (e.g., for GRB 940217) 03], deadtime effects and anomalous GRBs like GRB 941017 Our 
estimate uses the BATSE 4B fluence distribution, a Band function with Band j3 = —2.5 added to a > 30 MeV spectral 
component with a —2 photon index in the EGRET/GLAST range 1 12], and power-law approximations for the EGRET 
and GLAST effective areas fl3h . 

Fig.l, right panel, gives the resulting estimate of the number of GRBs GLAST should observe per year full-sky for 
a fluence ratio of 10%, for different integral photon numbers and energies. We find that there will be w 350 (« 25) 
GRBs/yr full-sky from which the GLAST LAT would detect > 5 photons with energy E > 30 MeV _[> 1 GeV), and 
« 25 (« 8) GRBs/yr full-sky with > 1 photon with E > 5 GeV (> 10 GeV). A similar conclusion Q, obtained by 
extrapolating BATSE results fitted with the Band function to LAT energies, gives an estimate of 50-70 GRBs/yr with 
> 5 photons with E > 30 MeV. Our analysis also shows that GRBs give very little ( ;$ 1% ) contribution to the diffuse 
extragalactic y-ray background. 
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